Epidemic spread in weighted scale-free networks 
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In this letter, we investigate the detailed epidemic spreading process in scale-free networks with 
links' weights that denote familiarity between two individuals and find that spreading velocity 
reaches a peak quickly then decays in a power-law form. Numerical study exhibits that the nodes 
with larger strength is preferential to be infected, but the hierarchical dynamics are not clearly 
found, which is different from the well-known result in unweighed network case. In addition, also 
by numerical study, we demonstrate that larger dispersion of weight of networks results in slower 
spreading, which indicates that epidemic spreads more quickly on unweighted scale-free networks 
than on weighted scale-free networks with the same condition. 

PACS numbers; 89.75.-k, 89.75.Hc, 87.23.Ge, 05.70.Ln 
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Many social, biological, and communication systems 
can be properly described as complex networks with ver- 
tices representing individuals or organizations and links 
mimicking the interactions among them. Recently, the 
ubiquity of a power-law degree distribution in real-life 
networks has attracted a lot of attentionP|. Examples 
of such networks (scale-free networks or SF networks 
for short) are numerous: these include the Internet, the 
World Wide Web, social networks of acquaintance or 
other relations between individuals, metabolic networks, 
integer networks, food webs, etc^]. The ultimate goal 
of the study on topological structure of networks is to 
understand and explain the workings of systems built 
upon those networks, for instance, to understand how 
the topology of the World Wide Web affects Web surf- 
ing and search engines, how the structure of social net- 
works affects the spread of diseases, information, rumors 
or other things, how the structure of a food web affects 
population dynamics, and so on. 

Recent studies on epidemic spreading in SF networks 
indicate a particular relevance in the case of networks 
characterized by complex topologies and very heteroge- 
neous structures that in many cases present us with 
new epidemic propagation scenarios ^^], such as absence 
of any epidemic threshold 0, hierarchical spread of epi- 
demic outbreaks!^, and so on. The new scenarios are 
of practical interest in computer virus diffusion and the 
spreading of diseases in heterogeneous populations. Fur- 
ther more, they also raise new questions on how to pro- 
tect the networks and find optimal strategies for the de- 
ployment of immunization resources 6]. However, so far, 
studies of epidemic spread just focus on unweighted SF 
networks, and a detailed inspection of epidemic spreading 
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FIG. 1: Density of infected individuals versus time in a BBV 
network with — 10*, 5 — 3.0, luq — 1.0 and m = 3, the four 
numerical curves i{t) correspond with parameter a =0.5, 0.4, 
0.33 and 0.25 respectively, form bottom to top. 



process in weighted SF networks is still missing while real 
networks, such as population and Internet, are obviously 
scale-free and with links' weights that denote familiarity 
between two individuals(like people or computers), re- 
spectively. One can easily take cognizance of how the 
links' weights affect the epidemic spreading process. For 
instance, if your little son gets flu, then you will be in- 
fected in all probability, since you two contact each other 
very frequently(i. e. of large familiarity) . By contraries, 
it is unlikely that you will be infected by your unfamiliar 
colleague just because of saying hello to him this morn- 
ing. 

In this letter, we intend to provide a first analysis of 
the time evolution of epidemic spreading in weighted 
SF networks. The weighted SF network model used 
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FIG. 2: Density of infected individuals versus time in a BBV 
network with A'^ = 10*, a — 2.0, ujo = 1-0 and m — 3, the four 
numerical curves i{t) correspond with parameter S =0.4, 0.3, 
0.2 and 0.1 respectively, form bottom to top. The inset shows 
the relationship between the dispersion of weight (fj.) and the 
value of S. 

in this letter is one of the most well-known model in- 
troduced by Barrat, Barthelemy, and Vespignani (BBV 
networks) Q , whose degree, strength and weight distribu- 
tions are power-law distributions with heavy tails. The 
BBV model suggests that two main ingredients of self- 
organization of a network in a weighted scale-free struc- 
ture are strength preferential attachment and weights' 
dynamics. These point to the facts that most networks 
continuously grow by the addition of new vertices, new 
vertices are preferentially attached to existing vertices 
with larger strength, and the creation of new links will 
introduce variations of the existing weight distribution. 
More precisely, the weight of each new edge is fixed to 
value uJo] if a new vertex linked to an existing vertex i, 
then the local rearrangement of weights between i and 
its neighbors j according to the simple rule 

^ LOij + AcJy (1) 

where 

AWy = (5^ (2) 
Si 

Si is the strength of node i, expressed by Si = J^j'^ij- 
This rule considers that the establishment of a new edge 
of weight uJo with the vertex i induces a total increase 
of traffic S that is proportionally distributed among 
the edges departing from the vertex according to their 
weights. Since BBV networks are of the same proper- 
ties (e.g. power-law distribution of degree, strength and 
weight) as many social networks (e.g. friendship net- 
works and scientists collaboration networks) and techni- 
cal networks (e.g. Internet and WWW), it is reasonable 
to investigate epidemic spreading on BBV networks. 

In order to study the dynamical evolution of epidemic 
spreading we shall focus on the susceptible-infected (SI) 



FIG. 3: Spreading velocity at each time f in a BBV net- 
work with TV = 10*, (5 = 3.0, uo ~ 1.0 and m = 3,when 
a = 0.33{square) and a — 0.25{circle). The inset shows 
the curves in log-log plot. The data are averaged over 200 
experiments. 

model in which individuals can be in two discrete states, 
either susceptible or infected 0. Each individual is rep- 
resented by a vertex of the network and the links are the 
connections between individuals along which the infec- 
tion may spread. The total population(the size of the 
network) N is assumed to be constant thus if S{t) and 
I(t) are the number of susceptible and infected individ- 
uals at time t, respectively, then N ^ S{t) + I{t) . In 
weighted networks, we define the infection transmission 
by the spreading rate, 

A,, = (^)",a>0 (3) 

at which susceptible individual i acquire the infection 
from the infected neighbor j, where a is a positive con- 
stant and ujm is the largest value of Wij in the network. 
Obviously, more familiar two individuals(i.e. with larger 
weight) may infect each other with greater probability. 
According to Eq.(3), one can quickly obtain the proba- 
bility that an susceptible individual i will be infected at 
the present time step is: 

A.(t)=i- n (1-^^^) (4) 

where Ni{t) is the set of all i's infected neighbors at time 
t. 

We start by selecting one vertex randomly and assume 
it is infected. The diseases or computer virus will spread 
in the networks in according with the rule of Eq.(3). In 
Fig.l, we plot the density of infected individuals versus 
time in a BBV network with N = 10'', 5 = 3.0, ljq = 1.0 
and m = 3. Since < 1, the smaller a is, the more 
quickly infection spreads. It is natural that larger value of 
5 induces larger dispersion of weight of networks. Then, 
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FIG. 4: Consider the cases of 5 < 1.0 and S — 1.0. (a) and (b) 
show the spreading velocity at each time f in a BBV network 
with iV = lO^.tJo = 1.0 and m = 3, when S = 0.5, a = 2.0 
and (5 = 1.0, a = 1.0, respectively. The inset shows the curves 
in log-log plot. The data are averaged over 200 experiments. 



a direct question is that how the value of 6 impacts 
epidemic spreading behavior. In Fig. 2, we show that 
epidemic spreads more quickly while 6 is smaller. In 
other words, larger dispersion of weight of networks re- 
sults in slower spreading. That means epidemic spreads 
more quickly on unweighted scale-free networks than on 
weighted scale-free networks with the same condition. 

Obviously, all the individuals will be infected in the 
limit of long time as limt_»oo*(0 = 1- For the sake of 
finding optimal strategies to protect individuals from be- 
ing infected, we will study the details of spreading veloc- 
ity at the outbreak moment. The spreading velocity is 
defined as: 



di{t) I{t) - I{t - 1) 



dt 



N 



(5) 



where i{t) = We account the number of newly in- 
fected vertices at each time step and report the spread- 



ing velocity in Fig. 3. Apparently, the spreading velocity 
goes up to a peak quickly that similar to the unweighted 
network cases 0, leaving us very short response time to 
develop control measures. Moreover, what's new and in- 
teresting, velocity decays following power-law form after 
the "peak time" . At the moment of infection outbreaks, 
the number of infected individuals is very small, as well as 
after a very long time from the outbreak, the number of 
susceptible individuals is very small. Thus when t is very 
small(close to zero) or very large, the spreading velocity 
is close to zero, one can see the corresponding simula- 
tion result in Fig. 3. One may think that the velocity 
follows power-law behavior just because of the extreme 
case of S > 1.0. Now we consider the case oi S < 1.0 and 
6 — 1.0. Fig. 4 shows spreading velocity at each time t in 
a BBV network with N — 10'', = 1-0 and m = 3, when 
6 — 0.5,0; — 2.0 and 5 = 1.0, a — 1.0, respectively. It is 
obvious that epidemic spreading behavior dose not show 
sensitive dependence on the parameter 6, the reason of 
that fact will be explored deeply in our future publica- 
tions. 

In order to give a more precise characterization of the 
epidemic diffusion through the weighted networks, we 
measure the average strength of newly infected vertices 
at time t, define as: 



I{t)~I{t-l) 



(6) 



where Is{t) is the number of infected vertices with 
strength s. Fig. ^ shows the average strength of newly 
infected vertices S^inf (0 as a function of time t, and the 
curves exhibit that -§inf(i) displays a power-law behav- 
ior for large Sintit) oc t~'^, which is remarkably dif- 
ferent from the clear hierarchical feature on unweighted 
networks 5]. 

It is explicit that the individuals with larger strength 
are much more dangerous when they are infected, rather 
than the ones with smaller strength, thus if one want 
to protect most individuals being infected, the suscep- 
tible individuals with larger strength must be protect 
foremost. In Fig. 5, one can find that the individuals 
with larger strength is preferential to be infected, which 
means there is little time leaving us to find the Large In- 
dividuals^^ and isolate them. Therefore, at the outbreak 
moment of disease or computer virus, the dense crowd or 
pivotal servers must be protected primarily. Of course, 
the outcome is not a good news for practical operators, 
but it may be relevant for the development of contain- 
ment strategies. 

In summary, we have studied epidemic spreading pro- 
cess in BBV networks, and the present results provide 
a clear picture of the infection propagation in weighted 
SF networks. The numerical studies show that spreading 
velocity Vini{t) and average strength of newly infected 
vertices S-mfit) present power-law time behavior for large 
t, which is remarkably different from infection propaga- 
tion in unweighted networks. Also by numerical study, 
we demonstrate that larger dispersion of weight of net- 
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FIG. 5: Behavior of average strength of the newly infected 
vertices at time t for SI model spreading in a BBV network 
with A'" = W^jS = 3.0, oJo = 1.0 and m — 3, the inset shows 
that Siniit) represents power-law behavior, Sinf(t) c>c t~'' . 

works results in slower spreading, which indicates that 
epidemic spreads more quickly on unweighted scale-free 
networks than on weighted scale-free networks with the 
same condition. These results indicate that not only 
the topological structures of networks but also the links' 
weights affect the epidemic spreading process. Further 
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more, the detailed study of behavior of average strength 
of the newly infected vertices may be relevant for the 
development of containment strategies. 

However, up to now, there are so many important and 
fundamental problems that puzzle us and haven't been 
referred to in the present letter. Some of them have been 
partially solved and will be publicized in further publi- 
cation, and others are still unanswered. At the end of 
this letter, we will list part of them. How to analyze 
the average density of infected individuals versus time 
at the outbreak moment in weighted SF networks, and 
how about the dynamic behavior after "peak time"? Is 
the mean-field theory appropriate to solve this problem? 
How to design a optimal containment strategy, and how 
about the effective for various strategies, such as to pro- 
tect vertices at random, to protect vertices purposefully, 
to cut off links at random, to cut off links purposefully, 
and so on? 
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